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INTRODUCTION 

The colossal magnetoresistance (OMR) of manganese oxides0 and their transport 
properties are usually studied in the framework of the double exchange (DE) Hamilto- 
nian or the more general ferromagnetic Kondo lattice model (KLM)n . The essence of 
these models is the high-spin configuration of e^-electron and t2g-coie electron spins due 
to a strong ferromagnetic Kondo exchange interaction K ~ leV. The kinetic energy 
in the partially filled Cg band is lowered when neighboring spins are aligned leading to 
a low-temperature ferromagnetic phase, while the high-T paramagnetic phase is disor- 
dered with a high resistivity. In these considerations the orbital degeneracy of the Mn 
Cg orbitals is usually neglected. Although this model provides an explanation of the 
phases necessary for a qualitative understanding of OMR, it has been stressed that the 
double exchange mechanism is not sufficient for a quantitative description, i 

Doped manganites are characterized by strong correlations and the complex inter- 
play of spin-, charge-, and orbital-degrees of freedom as well as the coupling to the 
lattice, e.g. via Jahn- Teller coupling, oi This complexity is directly evident from the 
large number of phases in a typical phase diagram. To quantify the different mecha- 
nisms it is helpful to analyze experiments where for certain parameters one or the other 
degree of freedom is frozen out. Important experiments in this respect are the very 
detailed investigations of the optical conductivity of Lai-ajSr^-MnOa by Okimoto et al.. 
10 These experiments (see Fig. ^ show (a) a pseudogap in (riuj) for temperatures 
above the Curie temperature (paramagnetic phase) with essentially linear in uo, 
and (b) the evolution of a broad incoherent distribution in the range < < 1.0 eV 
below Tc, which still grows at temperatures below Tc/10, where the magnetization is 
already close to saturation. Such a temperature dependence of (j{uj) over a wide energy 
region is quite unusual as compared with other strongly correlated electron systems 
near a metal-insulator phase boundary, i Interestingly there is in addition a narrow 
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Figure 1: Temperature dependence of the optical conductivity of Lai_a;Sr^Mn03 for 
X = 0.175 (reproduced from Okimoto et ai). A temperature independent background 
has been subtracted from the experimental data. 

Drude peak with width of about 0.02 eV and little weight superimposed to the broad 
incoherent spectrum. That is, the motion of charge carriers is dominated by incoher- 
ent processes but there is also a coherent channel leading to a narrow Drude peak. 
In view of the large energy scale of 1 eV it is plausible that the orbital degeneracy is 
the source of this incoherent motion, since the spin degrees of freedom are essentially 
frozen out. Our study of the degenerate Kondo lattice model will support this point of 
view. Further optical studies for various 3D manganites&S as well as for the bilayer 
system Lai.2Sri 8Mn2O70 confirm the presence of the large incoherent absorption in the 
ferromagnetic state. 

Although the importance of these experiments was recognized immediately, the few 
theoretical studiesi^HU were confined to simplified models or approximations, thereby 
ignoring important aspects of the full many body problem. 

The aim of this work is to show that in a model which accounts for the orbital 
degeneracy, yet assumes that the spins are fully polarized, the broad incoherent spectral 
distribution of (t{u), its increase with decreasing temperature, as well as the order of 
magnitude of cr(a;) at small u can be explained. Our calculation also accounts for a small 
and narrow Drude peak as observed by Okimoto et al. in the range uj < 0.05 eV. This 
is a clear indication of coherent motion of charge carriers with small spectral weight, 
i.e. in a model where due to the orbital degeneracy incoherent motion is dominant. 

We start our discussion with a generic Hamiltonian for the manganite systems, 
i.e. the ferromagnetic Kondo lattice model with degenerate e^-orbitals, and derive for 
the spin-polarized case an effective model which contains only the orbital degrees of 
freedom. This orbital model consists of a hopping term between the same and different 
orbitals a and f3 on neighbor sites and an orbital interaction. Renaming a = a where 
a =t or I the model maps on a generalized anisotropic t-J model. The usual t-J 
model known from the cuprates appears as a special case. Our derivation includes the 
3-site hopping processes, which appear as a natural consequence of the strong coupling 
limit. Although such terms do not influence the orbital order for integer filling, they 



are important for the proper description of transport properties in the doped systems. 

Because of the complexity of the orbital model we present here a numerical study of 
the orbital correlations and of the frequency dependent conductivity. Although there 
exist studies of the interplay of orbital and spin order at integer filling for LaMnOa, 
the effect of doping has not been considered so far. The finite temperature diag- 
onalization@'@ serves here as an unbiased tool to study the optical conductivity. The 
change of orbital correlations as function of doping and temperature is investigated by 
means of an optimized orbital basis (OOB). 

Furthermore we present results for the KLM without orbital degeneracy which clearly 
show that this model fails to explain the broad incoherent (t{uj) spectra observed in 
the ferromagnetic phase as T — 0. i'0 The KLM would instead lead to a sharp 
Drude absorption because of the ferromagnetic alignement, ^ i.e. to spinless fermion 
behavior in the low-temperature limit. This underlines the importance of the eg orbital 
degeneracy, although there are alternative proposals invoking strong electron-phonon 
coupling and lattice polaron formationSi'S to explain the incoherency of the absorption. 



MODELS 

We start from the generic ferromagnetic Kondo lattice model for the manganites, 
where the Mn e^-electrons are coupled to core spins Si (formed by the t2g orbitals) and 
a local repulsion U between the electrons in the two eg orbitals: 

H = Hhand + Hint + H Kondo- (l) 

Hband describes the hopping of e^-electrons between sites i and j, Hint the interaction 
between eg electrons and Hxondo the couphng of e^-spins and t2g-core spins. The eg 
electrons to have two-fold orbital degeneracy labelled by a Roman index (a,b) and two 
fold spin degeneracy labelled by a Greek index (a, cr'). Explicitly, 

H,and = E Etdljiaa + E (^J^L^J^^ + H.C.). (2) 
iacr (ij)afe(T 

The hopping matrix elements form a real symmetric matrix whose form depends on 
the choice of basis in orbital space and the direction of the i — j bond. For the present 
study we shall use \x) ~ — and z) ~ (3z^ — r^)/\/3 as basis for the eg orbitals. 
From the Slater-Koster rules foUowsiiil: 

^ "^^\ t^' --t(^ ^\ (3) 
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which allows for inter-orbital hopping in the xy-plane, where the upper (lower) sign 
distinguishes hopping along x and y direction. The hopping matrix elements are defined 
in terms of the double exchange t of d^z^-r^ orbitals along the c-axis. Here t = V^p^/Acdp 
is determined by the Mn-d 0-p hybridisation V^pa and the corresponding level splitting 
Aedp- All matrix elements are negative except t^^ along x-direction, which is positive. 
The level splitting AEi = E? — can be controlled by uniaxial pressure or in the case 
of layered compounds even with hydrostatic pressure, il We shall focus here on the 
orbital degenerate case, i.e. AEi = 0, however we will keep the orbital energy term for 
the derivation of the orbital model. 



The large local electron-electron repulsion U is responsible for the insulating behav- 
ior in the case of integer band-filling 



Hint = Uariia-iriiai + Uab ^ UiaUib + Jab ^ d\^d\^,„, dia^i diba ■ (4) 
in i icra' 

Here Ua, Uab and Jab are the intra- and inter-orbital Hubbard and exchange interactions, 
respectively. The relevant valences of Mn are Mn^"*" (S=2) and Mn"^"*" (S=3/2) as was 
pointed out already by Zener. i The Hamiltonian Hband + Hint was considered by 
Zaanen and Oles who showed that in general rather complex effective Hamiltonians 
result for transition metal ions with partially filled d shell near orbital degeneracy. 
Here a simplified approach is preferred, with the electrons of Mn ions forming core 
spins of S" = 3/2, and thus we restrict the electron interactions in Hint to the Cg bands. 

Even in the case with local orbital degeneracy there is one lower Hubbard band 
which is partially filled in the case of hole doping. With one e^-electron per site these 
systems are Mott-insulators. Although the Jahn- Teller splitting can lead to a gap in 
the absence of U, it is not the primary reason for the insulating behavior for integer 
filling, i 

The interaction of itinerant electrons with the (S=3/2) core spins is given by 

HRondo = -K V S^i ■ d\^^a^„idiaa' (5) 



leading to parallel alignment of the d-electron spin with the core-spin Si. Since the 
d-electron kinetic energy is favored by a parallel orientation of neighboring spins the 
ground state becomes ferromagnetic. □ 

In the low-temperature ferromagnetic phase we may introduce an effective Hamilto- 
nian which contains only the orbital degrees of freedom assuming a fully spin-polarized 
ferromagnet. The spin degrees of freedom can also be eliminated by high magnetic 
fields. This opens the possibility to investigate the orbital order independent of the 
spin degrees of freedom, and to shed light on the nontrivial question how the orbital 
order is changed upon doping with e^-electrons or holes. 

The resulting model has similarities to the t-J model, where the spin- indices a and 
a' are now orbital indices. To stress the difference we use for the orbital indices the 
letters a, 6 and Unitary transformation H^rb = e~^He^^ and the restriction to 

states without double occupancy leads to the orbital t-J model which has the following 
structure: 

Horb = E Kdldia + E (tpladib + H.c.) + H^- (6) 

with the constraint that each site can be occupied by at most one electron, i.e. d\a = 
dl^{l — fiia)- Here and in the following the index a denotes the orthogonal eg orbital 
with respect to orbital a. The orbital interaction if^^^ follows as a consequence of the 
elimation of doubly occupied sites with energy U r^Uab- Jab^ 

fj' = _i V V t'^'^ t'"' ( ^ I ^ ^ 

^ juu' aba/3 ^ + "'^j ~ "'^j+u ^ ^ "^j -^j+u' ^ 

[5/3,6 '^j+ua^j/g'^jfe'^j+u'a — 5/3,5 ^j+ua^jfe^j/^^j+u'a • (7) 

Here u, u' = (±a, 0) or (0, ±a) are lattice unit vectors. The orbital interaction H'^^f^ = 
H^^.^ + -f^oU consists of two types of contributions: (i) 2-site terms (u = u'), i.e. similar 
to the Heisenberg interaction in the standard t-J model, yet more complex because 



of the nonvanishing off-diagonal and (ii) 3-site hopping terms (u 7^ u') between 
second nearest neighbors. In the half-filled case, i.e. one e^-electron per site, only the 

2- site interaction is operative and may induce some orbital order. In the presence of 
hole doping, i.e. for less than one e^-electron per site, both the kinetic energy and the 

3- site contributions in H'^^f^ lead to propagation of the holes and to a frustration of the 
orbital order. 

The complexity of the model can be seen if we express the orbital interaction i^o^^ 
in terms of pseudospin operators T-^ = — riih), T-^ = dladib and T-~ = dlf^dia- where 
a =t and b =| denote an orthogonal orbital basis in the ec,-space. To be specific we 
shall assume here a{b) = z{x). Keeping only the two-site contributions, i.e. u = u', 
one obtains an anisotropic Heisenberg Hamiltonian for the orbital interactions: 

^:S = -|E [ ittt" + tf)i^i^i-T-Tp-t?f^i^iT^Tr+TrT^) 

(ij> 

+ i^' + ^)i^ini+TtT[)-t^^t^iT^T^ + TrTr) (8) 
- (t?,«t?,^ - t'^h'^^) {T.'T.+ + T.'Tr ) - - t^h?,^) (T+r^ + T.-T.^) . 

J o o o o o J o J o o o 

The orbital interaction is equivalent to that studied by Ishihara et al. iB Since 
we are interested here in the effect of holes on the orbital order we will generally base 
our study on H'^^f^ (0) which includes the 3-site hopping processes as well. 

We note that in the special case t^'' = 6abt and E?- = equations (@) and (|^) are 
identical to the usual t-J model ( with J = At^ /U and including 3-site processes). For 
the orbital model we shall adopt the same convention for the orbital exchange coupling, 
i.e. J = AtyU. 

The orbital degrees of freedom in combination with strong correlations (i.e. no 
doubly occupied sites are allowed) is expected to lead to incoherent motion of holes, 
although quasiparticle formation with reduced spectral weight is a plausible expectation 
on the basis of what is known about the usual t-J model. 

In this paper we focus on the orbital degenerate case with Ef = E\ = 0. A detailed 
study of the influence of a finite level splitting will be presented elsewhere. 



OPTICAL CONDUCTIVITY 

We study the charge transport by calculating the optical conductivity 

aoi^uj) = 2TTe^D^5{uj) + a{uj). (9) 

The frequency dependent conductivity consists of two parts, the regular finite frequency 
absorption crUiA and the 5-function contribution which is proportional to the charge 
stiffness Z}J5§E3 . The latter vanishes in insulators. This contribution is broadened into 
a usual Drude peak in the presence of other scattering processes like impurities which 
are not contained in the present model. The finite frequency absorption (or regular 
part) cr{u) is determined by the current-current correlation function using the Kubo 
formula 

a(^) = — Re dte'-\Ut)j,). (10) 

JMuj Jo 

For the derivation of the current operator we introduce twisted boundary conditions 
via Peierls construction 

tg,j(A) = tS„jexp(-^| j^'^"i(x)crx). (11) 



Assuming A = {A^, Ay, 0) we obtain from the kinetic energy operator the x-component 
of the current operator = dH{Ax) /dAx. 



j+uab 



An additional contribution follows from the 3-site term in Eq.(6) 
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where the expression 0^^(j,u, u') is an abbreviation for the operator within the last 
brackets of Eq.(6). 

We stress here that the orbital order in the undoped case is determined exclusively by 
the (two-site) orbital interaction (7), while the 3-site processes in (6) only contribute in 
the doped case. Nevertheless the motion of holes and transport in general is influenced 
in a significant way by the 3-site current operator In fact from studies of the t-J 
model it is known that cr(co') is not only quantitatively but even qualitatively changed 
by these terms. 00 Although the 3-site terms are usually ignored in studies of the 
t-J model they are an important part of the strong coupling model and should not be 
droped in studies of the conductivity. One of the aims in our study of the orbital model 
is to analyse the effect of these 3-site terms in the model. 

The charge stiffness Dc can be determined in two ways: (a) using Kohn's relation, 
or (b) via the optical sum rule which relates the integrated spectral weight of the 
real part of the conductivity to the average kinetic energy: 



Here H'^^^ contains two contributions: (a) the usual kinetic energy ~ t (5) and (b) 
the 3-site hopping processes ~ t'^/U in Eq. (6). The case with off-diagonal hopping 
considered here is a generalization of the sum rules for the t-J modefi and for the t-J 
model including 3-site terms. Together with (8) this implies 



In the following we shall use Dc = S — Suj a.s abbreviation for this equation, where S 
denotes the sum rule expression and the finite frequency absorption. 

FINITE TEMPERATURE LANCZOS METHOD (FTLM) 

For the calculation of cr{uj) we use a generalization of the exact diagonalization tech- 
nique for finite temperature developed by Jaklic and Prelovsek. c3 In this approach the 
trace of the thermodynamic expectation value is performed by a Monte-Carlo sampling. 
The current-current correlation function in (9): 




(14) 




(15) 



C{uj)=Re dte^\Ut)jx). 
Jo 



(16) 



can be rewritten by introducing a complete set of basis functions | r) for the trace, and 
eigenfunctions | and | of if with eigenvalues and Ej (here r is an irrelevant 
label, which will get its meaning below): 

Ciu;) ^ ^ E E (r I I J. I ^ ^(^ • I J. I + El - i?}) (17) 

r=l i,j=l 

with 

^-EEe-^'^'^Krl^Dr- (18) 

r=l i=l 

These expressions are exact if one uses complete sets. The FTLM is based on two 
approximations when evaluating Eqs. ([T7|) for C{uj) and ([T8| ) for the partition function 
Z: (1) The trace is performed over a restricted number R of random states | r), and (2) 
the required set of eigenfunctions of H is generated by the Lanczos algorithm starting 
from the initial states | $q) =| r) and | $q) = jx \ r) / ^ {r \ \ r), respectively. The 
latter procedure is truncated after M steps and yields the relevant intermediate states 
for the evaluation of C^u). 

Detailed tests0 have shown that the results get very accurate already if i?, M <C iVg, 
where Nq is the dimension of the Hilbert space. The computational effort is significantly 
reduced because only matrices of dimension M x M have to be diagonalized, i.e. much 
smaller than the dimension Nq x Nq of the full Hamiltonian H. In practice typical 
dimensions are R = 100 and M = 100. 

A complication in the case of the orbital t-J model is that T^^^ does not commute 
with the Hamiltonian Horb (iD- Therefore the calculations for the orbital model are 
restricted to smaller clusters than in the t-J case, where different S^^^ subspaces can be 
treated separately. 



RESULTS 

Kondo-Lattice Model without Orbital Degeneracy 

We begin our discussion of the frequency dependent conductivity with results for 
the single orbital Kondo-lattice model excluding double occupancy 

-f^ = - E ^ii^Ldja - Jn'^Si- dl^aaa'dia'- (19) 

ijcr iaa' 

Nearest neighbor hopping = t is assumed and the actual calculations were performed 
for 5* = 1 core spins for different doping concentrations x. 

In Fig.0 the temperature dependence of the optical conductivity is shown for a 
chain at doping x = 2/6. The conductivity is given in dimensionless form apQ, where 
Pq = ha/e^ sets the dimension of the resistivity in a 3D s.c. crystal structure with lattice 
constant a. At strong coupling Jh ^ t the states of the model consist of a low-energy 
band corresponding to the S* + 1/2 configurations, and of the exchange-splitted bands 
separated by an energy ^Jh- The intraband transitions are shown in the large window. 
They give rise to an incoherent low-energy absorption extending up to the energy of 
the free carrier band-width ~ At. The spectrum of the high to low spin transitions is 
centered around uj ^ \ Jh and is shown in the inset. 

At low temperatures, however, when ferromagnetic correlations extend over the 
whole system, the low-energy spectrum is dominated by a narrow coherent peak at 
cij = 0, and the interband transitions vanish. At T = all weight is in the Drude peak 
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Figure 2: Frequency dependent conductivity <j{uj) of the single orbital KLM ior a six 
site chain with two electrons and Jh = Qt for different temperatures. Inset shows the 
spectrum of excitations to the exchange-sphtted band. Note that for T — > all spectral 
weight is in the Drude peak (From Jaklic, Horsch and MackEl). 
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Figure 3: Temperature dependence of the kinetic energy and of the optical spectral 
weights of the ferromagnetic Kondo-lattice model in the absence of orbital degeneracy 
for a 6 site chain with two holes. Here S = {H^l^) /2Nt and = J a{uj)duj/7ieH. The 
spectral weight of the interband optical transitions into low-spin states is measured by 
S2 = J a2{ijj)duj /neH. 




Figure 4: Optical conductivity cr{uj) calculated for the orbital t-J model for a N=10 
site plane with A^e=8 electrons and J=0.25t (U=16t) at different temperatures. The 
calculation includes the 3-site terms. The inset shows the optical conductivity cro{uj) 
including the Drude-part. The spectra were broadened with r=0.1t. 

(see Fig.|^). As the temperature is increased the spectral weight is transfered from the 
coherent part to the broad incoherent spectrum and to the exchange-split excitations. 
This clearly shows that the model without orbital degeneracy does not explain the 
anomalous absorption in the ferromagnetic state of manganites (Fig.|I]). 

Optical Conductivity for the Orbital t-J Model 

Typical results for the frequency and temperature dependence of cr{uj) for a two- 
dimensional 10-site cluster with 2 holes, i.e. corresponding to a doping concentration 
X = 0.2, are shown in Fig. ^ with and Fig. ^ without 3-site hopping terms, respec- 
tively. The cr{uj) spectra are bell-shaped and increase with decreasing temperature. The 
width of the distribution measured at half-maximum is uJi/2 ~ 2.5t and is essentially 
independent of temperature. Calculations for a 3D cluster yield a rather similar cr{uj) 
distribution with a slightly larger value for the width uJi/2 ~ 3.5t. 

If we compare this latter (3D) value with cu^Jg ~ 0.7 eV found from the data of 
Okimoto et al. (for x = 0.175 and T = 9K), we obtain an experimental estimate for 
the parameter t: t"^* ~ 0.2 eV. For comparison a rough theoretical estimate based on 
Harrison's solid state table yields t = V^p^/Ae ~ 0.4 eV. Hence the orbital t-J model 
explains in a natural way the ~ 1 eV energy scale of the incoherent absorption in the 
experiments. Although shows the anomalous increase with decreasing tempera- 
ture, we do not claim here that the orbital model accounts for the full temperature 
dependence in the ferromagnetic phase (T < Tc). The description of the complete 
T-dependence certainly requires the analysis of the degenerate-orbital KLM, i.e. the 
additional inclusion of the spin degrees of freedom. Yet we believe that the results can 
be compared with the data by Okimoto et al. in the low-T limit, i.e. in the regime 
where the magnetization is saturated. 




Figure 5: Optical conductivity as in Fig.^ but without 3-site contributions. Due 
to the stronger incoherence there is a smaller Drude peak (inset) and a smaller DC- 
conductivity. (From Horsch, Jaklic and MackS). 



To compare with experimental data we present the conductivity in dimensionless 
form and with the proper dimensions of a 3D conductivity, assuming that the 3D system 
is formed by noninteracting layers. For a cubic lattice po = ^a/e^. If we consider 
Lai-xSr^MnOa with lattice constant a = 5.5 A and h/e^ = A.llKQ we obtain po = 
0.23 ■ 10"'^ Qcm. Hence the low frequency limit of the regular part of the conductivity 
(j(a; 0) = o"o"^ ~ 0.3—0.4-10^ (ficm)"^. This is consistent with the order of magnitude 
for the low frequency limit of the incoherent part of the cr(a;) data of Okimoto et al. 
for Lai-xSr^^MnOs cx™'^ ~ 0.4 ■ 10^ (0.3 ■ 10^) (ficm)"^ for the doping concentrations 
X = 0.175 (0.3), respectively, at T = lOK. 

Hence we conclude that besides the energy scale also the absolute value of the in- 
coherent part of the experimental cr{u}) spectrum is consistent with the orbital model. 
We stress that the value (order of magnitude) of a^'^ is essentially fixed by the con- 
ductivity sum rule and the scale uJi/2, as long as the conductivity is predominantly 
incoherent. The insets in Fig.^ and ^ show cro(co'), Eq.(^, with the Drude absorption 
at low frequency included, where we used an ad hoc chosen parameter F = O.lt to 
broaden the 5-function. This value corresponds to the experimental width F ~ 0.02 
eV taken from the (small) Drude peak observed by Okimoto et al. i While our model 
yields the weight of the Drude peak D^, it does not give F, which is due to extrinsic 
processes (e.g. scattering from impurities and grain boundaries). The charge stiffness 
Dc together with F determines the DC- conductivity. 

Adopting the experimental estimate for F we find for the example in Fig.^ (inset) 
for the low temperature DC-conductivity aoc ~ 6.5-10^ (f2cm)~^ and for the resistivity 
p ~ 0.15-10"^ Qcm, respectively. For comparison, the experimental range of resistivities 
is e.g. 0.1 — 1.0 ■ 10^^ Qcm in the ferromagnetic metallic phase of Lai-^^Sr^-MnOaiil . 
As we shall see below, the Drude weight and therefore the DC- conductivity depend 
considerably on the value of the exchange coupling J, and whether the 3-site hopping 
processes in the model are taken into account or not. These terms have a strong effect on 
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Figure 6: Temperature dependence of the kinetic energy S (solid line), the incoherent 
spectral weight (dotted) and the Drude weights Dc (dash-dotted line) and DKohn 
(diamond) for a N=10 site plane with Ne=8 electrons and J=0.25t (U=16t) in units of 

the coherent motion of charge carriers. After this discussion it should be clear, however, 
that the DC-conductivity at low temperatures is largely determined by extrinsic effects. 

Figure ^ shows the temperature variation of the sum rule (kinetic energy), the finite 
frequency absorption 5*^ and the Drude weight for the model including the 3-site 
terms. Whereas at high temperature the sum rule is essentially exhausted by the finite 
frequency absorption (t{uj), we find at low temperatures a significant increase of the 
Drude weight. In this particular case Dc contributes about 40% to the sum rule at 
T/t = 0.3. 

The conductivity data shown in Figs. ^ and ^ for T/t > 0.5 is characteristic for 
a 2D orbital liquid state, which is stabilized by thermal fluctuations. The significant 
increase of coherency in Fig. |^ below T/t = 0.3 is due to the onset of x^-y^ orbital order 
in the planar model (see discussion below). This ordering is characteristic for the 2D 
version of the orbital t-J model and does probably not occur in cubic systems for which 
an orbital liquid ground state was proposed. Therefore we consider the T/t = 0.3 
data for Dc to be more appropriate for a comparison with the low-temperature 3D data 
than the T = values, which are enhanced due to orbital order. 

We have also calculated Dc for T = using Kohn's relationil 
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(20) 



which yields a consistent value (Fig.||). Here $ = eaAa/h is the Peierls phase induced by 
the applied vector potential with component A^,. For the evaluation of Eq.|2y we followed 
Ref.t3 and first searched for the minimum of the ground state energy Eq with respect to 
the vector potential, 123 and then calculated the curvature. Since this procedure is quite 
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Figure 7: Kinetic energy S versus J for a N=10 site 2D cluster with N(,=7 and 8 elec- 
trons, with (full) and without (open symbols) the 3-site terms {T /t = 0.3), respectively. 
The dash-dotted line (circles) represents a 16 site cluster with iVe=14 and without 3-site 
terms (From Horsch, Jaklic and MackEj). 



cumbersome, we have determined most of our data via the sum rule. The resulting 
sets of data show systematic trends as function of J and doping. 

From these results we expect in analogy with the t-J model that the single-particle 
electron Green's function is characterized by a pronounced quasiparticle peak to explain 
such a large fraction of coherent transport. 

From Fig.|^ we also see that the weight under a{uj), i.e. S^, does not further increase 
for temperatures below T = 0.5t, although the temperature variation in Fig.^ seems to 
suggest a significant further increase towards lower temperatures. Conductivity data 
for T < 0.5t is not shown in Figs. ^ and ^ because it shows pronounced discrete level 
structure, and probably requires larger clusters for a careful study. Nevertheless it 
appears that (t{u!) develops a pseudogap in the orbital ordered phase. H Integrated 
quantities on the other hand are much less influenced by such effects. 



The Role of 3-Site Terms 



In the following we wish to shed more light on the role of the 3-site processes in the 
orbital Hamiltonian and in the current operator. The importance of the 3-site term 
becomes particularly clear from the J-dependence of the sum rule S (Fig.|^) and the 
charge stiffness Dc (Fig.|^) taken at temperature T/t = 0.3. These low temperature 
values of S and D^. are approximatively independent of the exchange parameter J in 
the absence of 3-site terms. A similar observation was made for the t-J model, illii 
Moreover one can see that the sum rule S is proportional to the doping concentration x 
for the 3 cases shown (x=0.125, 0.2, and 0.3). When 3-site terms are taken into account 
both S and acquire a component which increases linearly with J. These general 
features are fully consistent with results for the t-J model. 11^ 
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Figure 8: Drude weight versus J calculated via the sum-rule (|T5|) for a N=10 
site planar cluster with Ne=7 and 8 Cg electrons, with and without the 3-site terms, 
respectively. 



The change due to the 3-site terms is particularly large for the charge stiffness, 
which defines the weight of the low frequency Drude peak. The relative weight in the 
Drude peak Dc/S is shown in Fig. | as function of J = At'^/U. Small changes in J lead 
to considerable changes in the Drude weight and the coherent motion of carriers. For 
J = 0.25 (0.5) the 3-site terms lead to an increase of Dc by a factor of 3 (5), respectively, 
for the 8 electron case. 

When we take the value t ~ 0.2 eV, which we determined from the comparison 
of a{uj) with the experimental data, and [/ ~ 3 eV, we obtain U/t ~ 15 and 
J/t = At/U ~ 0.25 (J = 0.05 eV). For such a value for J we expect a relative Drude 
weight Dc/S ~ 0.4, which is larger than the experimental ratio D^/Neff ~ 0.2 found 
by Okimoto et al. for x = 0.175, identifying here the effective number of carriers iVe//0 
with 5*. The small value for Dc/N^ff found in the experiments suggests that J is not 
larger than the value estimated, otherwise we would expect a too large relative Drude 
weight. A larger value for U would also have the effect to reduce Dc- However, we have 
to stress here that certainly calculations on larger clusters must be performed, before 
possible finite size effects in can be quantified. 



Comparison with the standard t-J Model 



It is interesting to compare these results for the orbital t-J model with those ob- 
tained for the usual t-J model describing the carrier motion in the copper-oxygen planes 
in high-Tc superconductors. We note that the t-J model is a special case of the orbital 
model with taa = tbb and tab = 0- Figures | and show results for a 4 x 4 cluster with 
Nh = 3 holes. The results for (j[u!) show for both models a large incoherent absorption, 
which is bell-shaped in the case of the orbital model, while in the t-J case there is a 
continuous increase of <j{uj) as u approaches zero. The temperature dependence of the 
sum rule, the incoherent and the Drude weight shown in Fig. [1^ are quite similar for 
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Figure 9: Optical conductivity for the standard t-J model (taa = thh = i and tab = 0) at 
temperatures T = 0.5, 1.0, 2.0 and 4.0 which should be compared with the corresponding 
data for the orbital model. Here Nh/N = 3/16 and J/t = 0.4. 

both models, with the important exception that the orbital model goes into the ordered 
phase at low temperatures (T < 0.3) which is accompanied by a marked increase of 5* 
and Dc (Fig.|). A more detailed study of the T-dependence may be found in the review 
by Jaklic and Prelovsek. @ 

Comparison with other Work 

The optical conductivity for the orbital non-degenerate Kondo lattice model was 
studied by Furukawa using the dynamical mean-field approximation. An important 
result of this calculation is, that the weight of intraband excitations within the lower 
exchange-split band should be proportional to the normalized ferromagnetic magnetiza- 
tion M/Msat- Experimentally, however, Neff{^ S) is still increasing when temperature 
is lowered even though M is already saturated. This contradiction to the prediction of 
the simple double-exchange model implies that some other large-energv-scale scattering 
mechanism survives at low temperature, where the spins are frozen. The dynamical 
mean field theory yields one low energy scale, and not two, i.e. there is no Drude peak 
plus incoherent structure. 

Shiba and coworkers et a/.0 approached the problem from the band picture. They 
analysed the noninteracting two-band model (2) (i.e. without constraint) and argued 
that the interband transitions within the Cg orbitals may explain the anomalous absorp- 
tion in the ferromagnetic phase of Lai_2,Sr2,Mn03. While the frequency range of the 
interband transition is found consistent with the anomalous absorption, the structure 
of cr{u) differs. In particular the noninteracting model leads to cr(u;) ~ u for small cu. 
Moreover they calculated the ratio S^/D^ ~ 0.85, i.e. there is more weight in the Drude 
peak than in the regular part of the conductivity a{uj). This ratio is found to be rather 
insensitive to doping in the range 0.175 < x < 0.3. This result differs considerably from 
the experimental data of Okimoto et al. where this ratio is about 4 for x = 0.175. 
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Figure 10: Temperature dependence of the kinetic energy S (solid line), the incoherent 
spectral weight S^^ (dotted) and the Drude weight Dc (dash-dotted line) for the t-J 
model. The data was calculated for a 2D system with 16 sites and 3 holes for J/t=OA. 

Orbital fluctuations in the ferromagnetic state and their effect on the optical conduc- 
tivity were studied by Ishihara et alS^ using slave fermions coupled to bosonic orbital 
excitations. Ishihara et al. deal with the concentration regime xt ^ J, where the 
kinetic energy is expected to dominate the orbital exchange energy, and arrive at the 
conclusion that the quasi two-dimensional nature of the orbital fluctuations leads to an 
orbital liquid in (3D) cubic systems. The orbital disorder is treated in a static approxi- 
mation in their work. The cr(a;) spectrum obtained by this approach for T = O.lt has a 
quite similar shape as the experimental spectrum at low temperature (with ujii2 ~ 3t), 
however, because of the assumed static disorder, there is no Drude component. A 
further problem with the orbital liquid, as noted by the authors, is the large entropy 
expected for the orbital disordered state which seems to be in conflict with specific heat 
measurements. 

More recently a microscopic theory of the optical conductivity based on a slave 
boson parametrization which combines strong correlations and orbital degeneracy was 
developed by Kilian and Khaliullin§ for the ferromagnetic state at zero temperature. 
Their approach yields a highly incoherent spectrum due to the scattering of charge 
carriers from dynamical orbital fluctuations, and a Drude peak with strongly reduced 
weight. However (t{uj) is depressed at low frequency in this calculation for the orbital 
model and an additional electron-phonon mechanism is invoked to obtain results similar 
to the experimental spectral distribution. This theory further accounts for the small 
values of the specific heat and therefore supports the orbital liquid scenario for the 
ferromagnetic cubic systems, lil 

We would like to stress that the proposed orbital liquid statelll is a property of cubic 
systems, and not a property of the 2D or quasi-2D versions of the orbital model. In 
the planar model the cubic symmetry is broken from the outset and holes presumably 
cannot restore it. This is different from the physics of the t-J model, where holes restore 
the spin rotational symmetry (which is spontaneously broken in the antiferromagnetic 



ordered phase), because this symmetry is respected by the t-J Hamiltonian in any 
dimension. 

For a deeper understanding of the results for the optical conductivity, we shall 
analyse below the structure of the orbital correlations in the doped and undoped two- 
dimensional orbital model in detail. 



DOPING DEPENDENCE OF ORBITAL CORRELATIONS 

In the absence of holes the interaction @) between orbital pseudo-spins will lead 
to an orbital order below a certain temperature ~ J due to the anisotropy of the 
Hamiltonian. Doping will destroy this order and lead either to a disordered orbital 
liquid, or may generate a new kind of ordered state which optimizes the kinetic energy 
of the holes. To investigate this question for an anisotropic model it is in general not 
sufficient to calculate simply e.g. the correlation function < T-^T? > defined with respect 
to the original orbital basis, because the relevant occupied and unoccupied orbitals 
may be different from the original orbital basis chosen. In the following we therefore 
introduce a local orthogonal transformation of the orbitals on different sublattices by 
angles and ■?/', respectively. On the A-sublattice: 

\z) = cos{(f))\z) + sin{(j))\x) 

\x) = —sin{(j))\z) + cos{(f))\x), (21) 

and a similar rotation with angle ip on the B-sublattice. This amounts to new operators, 
e.g. 

= cos{2<P)Tl + sin{2ct>)T: (22) 
f^^^ = cosmT^_,^ + szn{2ij)Tt^^ (23) 

and new correlation functions, e.g. (T^T^) is defined as 

{f^fl^j,) = cos{2cl>)cos{2i,){TlT:^^) + sin{2cl>)sin{2i^){T^T^^^) 

+cos{2ct^)sin{2^){T^T:^^) + .m(20)cos(2^)(7;-7;%) (24) 

The angles and ip are now chosen such that the nearest neighbor correlation function 
(T-^T-^) takes a maximal (ferromagnetic) value. With this convention the orbital order 
is specified by the corresponding local quantization axis, that is by the values of the 
rotation angles and ip. We call this the optimized orbital basis (GOB), ill 

In the undoped 2D case (x-y plane) these angles are (p = 45° and ip = 135° cor- 
responding to an :^(|a;) + \z)) and -^Hx) — \z)) order at low temperatures (see Fig. 

|lT] and |l^). (T^Tl^j^) takes a value only slightly smaller than 0.25 at low tempera- 
ture (Fig. ^), which implies that quantum fiuctuations are small in this state. Hence 
the orbital correlation function (T/Tj^^) of the undoped orbital model is alternating 
(i.e. antiferromagnetic in the pseudo-spin language for the orbital degrees of freedom). 
That the AF-order seems to be established at a finite temperature for the small clusters 
means that the correlation length ^ gets larger than the system size L. At small tem- 
peratures the correlation functions have about the same value independent of R, i.e. 
the data shows no spatial decay. Above T/t = 0.1 the orbital correlations are small and 
show a pronounced spatial decay, which can be considered as a signature of a thermally 
disordered orbital liquid state. 
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Figure 11: Contourplot of the rotated nearest neighbor orbital correlation function 
{TfTl^^), R={1,0), as function of the angles (p and ip for a N=10 site planar cluster 
with A^e=10 electrons, J=0.25t (U=16t) and T=0.01t. White regions correspond to 
positive (ferromagnetic) and black areas to negative (antiferromagnetic) correlation 
functions, i.e. {T/T^^^)>0.23 (<-0.23), respectively. The 25 contour lines are chosen 
equidistant in the intervall [-1/4, +1/4]. (From Horsch, Jaklic and Mack0) 



Experimentally the coexistence of antiferromagnetic (or staggered) orbital order anc 
ferromagnetic spin order has been established in the low-doping regime of LaMnOa. 
It should be recalled that LaMnOs has an A-type antiferromagnetic spin structure, ^ 
where ferromagnetic layers are coupled antiferromagnetically along the c-axis. In the 
orbital ordered case the MnOg octahedra are deformed, and Mn Sd^^'^^r^ and 3d3y2_j.2 
orbitals with orientation alone the x- and y-axis, respectively, have been considered as 
relevant occupied orbitals. e3e3 For the planar model we find here the alternate occu- 
pation of ■:^{\x) — \z)) and + \z)) orbitals (Fig. |12[) which are also oriented along 
the X- and y-axis, however they differ with respect to their shape (and have pronounced 
lobes along z). 

Doping of two holes in a ten site cluster is sufficient to remove the alternating orbital 
order and establishes a ferromagnetic orbital order {(p = ip = 0) as can be seen from 
Figs. and 0. This corresponds to preferential occupation of orbitals with symmetry 
x"^ — as shown in Fig.^. The pseudospin alignment can be considered as a kind of 
double exchange mechanism in the orbital sector. An interesting feature in Fig.[l^for the 
doped case is the fact that the correlations do not show significant spatial decay even at 
higher temperatures where the correlations are small. Moreover the correlations in the 
doped case appear to be more robust against thermal fiuctuations than in the undoped 
case. The characteristic temperature (determined from the half-width in Figs. 13 and 
rBj) is T* ~ 0.2t for the 2-hole case while for the undoped system T* ~ 0.05t. This 



trend is consistent with the fact that in the doped case order is induced via the kinetic 
energy and the corresponding scale t is larger than J = 0.25t. Although a more detailed 
analysis would be necessary to account properly for the doping dependence. 

It has been shown recently by Ishihara et a/.il that in the layered manganite com- 
pounds hydrostatic pressure leads to a stabilization of the x'^-y'^ orbitals as well. 
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Figure 13: Temperature dependence of (T/Tj!,_^) with cj) = 45° and ip = 135° for 
nearest-neighbors and ip = 45° for next-nearest neighbors, respectively. Results are 
shown for different distances R for an undoped N = 10 site planar cluster. Parameters 
as in Fig.|lT]. The strong increase for T/t < 0.1 is due to the onset of alternating orbital 
order. (From Horsch, Jaklic and MackEll) 
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Figure 14: Contourplot of {T^T^j^j^ for R = (1,0) as function of angles and ip for a 
N=10 site cluster with two holes. Otherwise same parameters as for Fig.pT. 




Figure 15: x'^-y'^ orbital order in the doped phase. 
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Figure 16: Temperature dependence of (T^Tl^j^) with = 0° and ip = 0" for a 10-site 
cluster with two eg holes. Parameters as in Fig. 11. The strong increase below T/t = 0.3 
signals the onset of x^-y^ orbital order in the planar model. 



A comparison of the doping dependence of spatial correlations with the t-J model is 
given in Fig.|l^ In the undoped case the orbital model shows almost classical (orbital) 
Neel order, whereas in the t-J model correlations are strongly reduced by quantum 
fluctuations. The T = correlation function of the t-J model shows long-range order 
(LRO) consistent with a strongly reduced sublattice magnetization ~ 0.3 (in the 
thermodynamic limit). 

In the orbital model the effect of the two holes {x = 0.125 and 0.2) is strong enough 
to induce ferromagnetic (orbital) correlations, i.e. with a prefered occupation of x^-y^ 
orbitals. In the t-J model instead spin-correlations decay rapidly for two holes on a 
20-site cluster (x = 0.1), which is consistent with the notion of an AF spin-liquid state 
(see inset Fig.p!7|). 

The preference of ferromagnetic orbital order in the doped case is due to the fact 
that t"^ depends on the orbital orientation, and t"'"' 3> t'^^. One would expect that for 
sufficiently strong orbital exchange interaction J = At^ /U one reaches a point where 
antiferromagnetic orbital interactions and ferromagnetic correlations due to the kinetic 
energy are in balance. This quantum critical point between FM and AF orbital order 
turns out to be at quite large J values. For the two-hole = 0.2) we have found 

that this crossover happens for quite large orbital interaction J ~ 2, i.e. at a value 
about an order of magnitude larger than typical J values in manganite systems. 

Before closing this section, we shall explain the different orbital occupancy in the 
doped and undoped case in more physical terms. The overlap of the two sets of orthogo- 
nal Cg orbitals on neighboring sites changes with the angles and if). (1) In the undoped 
state the direct hopping between the (predominantely) occupied orbitals -^{\x) + \z)) 
and — \z)) on neighboring sites is small (Fig.|T2|). The hopping between these 

occupied orbitals is blocked anyhow because of Pauli's principle. However t"-^ between 
an occupied orbital on one site and an unoccupied orbital on a neighbor site is extremal. 
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Figure 17: Orbital correlation functions for T = as function of distance (squared) for 
a 10-site and a 4 x 4 cluster with zero and two holes {J/t = 0.25). In the undoped case 
the correlation function is antiferromagnetic, whereas in the doped case orbital correla- 
tions are ferromagnetic. Antiferromagnetic orbital order is indicated by full (negative) 
and open (positive) symbols, respectively. Inset: The spin correlations (S'f S'f_,_j:j) for a 
20-site 2D Heisenberg model and the t-J model with two holes are shown for compari- 
son [t-J data for J = 0.4 taken from Horsch and Stephan0). The two-hole case shows 
a rapid decay of correlations characteristic for an antiferromagnetic spin liquid. 



This leads to a maximal antiferromagnetic orbital exchange interaction, and thereby to 
a lowering of the energy. (2) In the doped case instead, the x^-y"^ orbital occupancy is 
prefered on all sites (Fig.|15]), which implies a large hopping amplitude for holes in the 
partially occupied x^-y^ band, whereas the off-diagonal hopping matrix element t"* and 
the orbital interactions are reduced for this choice of occupied orbitals. Hence the main 
energy gain in this case is due to the correlated motion of holes in the x^-y^-band. Nev- 
ertheless the off-diagonal hopping plays an important role in the x^-y^ ordered state, 
and leads to the strong incoherent features characteristic for the conductivity a{uj). 

Finally we note that recent experiments by Akimoto et a/.S proved the existence 
of an A-type antiferromagnetic metallic ground state in a wide doping concentration 
regime of the 3D system (Lai_2Nd^)i_a;Sra;Mn03 and in particular for Lai-a^Sr^^Mn O3 
at high doping x = 0.52 — 0.58. The latter compound was investigated by Okimoto 
et al. for smaller doping x where the ground state is the uniform ferromagnetic state. 
The A-type antiferromagnet consists of ferromagnetic layers which are coupled anti- 
ferromagnetically. Akimoto et al. propose that x^ — y"^ orbitals should be occupied 
in this state and form a highly anisotropic 2D band. This seems consistent with our 
finding of — orbital order in the ferromagnetic planar system where the cubic 
symmetry is explicitely broken. An experimental study of the optical properties in this 
concentration regime at low temperature would be particularly interesting. 

SUMMARY AND CONCLUSIONS 

Our work was motivated by the experimental study of the frequency dependent 
conductivity of Lai_a;Sra;Mn03 by Tokura's group, i In particular these experiments 
show in the saturated low-T ferromagnetic phase a broad incoherent spectrum and at 
small frequency in addition a narrow Drude peak, which contains only about 20% of the 
total spectral weight. This experiment by itself shows that the ferromagnetic metallic 
state (although fully spin polarized) is highly anomalous . These features cannot be 
explained by the standard double exchange model, but require an additional mechanism. 
The most natural mechanism, one can think of, is the twofold degeneracy of the Cg 
orbitals. 

We have studied the optical conductivity, the charge stiffness and the optical sum 
rule for an effective Hamiltonian that contains only the e^-orbital degrees of freedom. 
This model follows from the more general Kondo lattice model, when the spins are 
ferromagnetically aligned. The 'orbital t-J model' is expected to describe the low- 
temperature physics of manganese oxides in the ferromagnetic saturated phase. 

Our results for cr(a;) in the 2D orbital disordered regime show that these character- 
istic features follow from the orbital model. We believe that the optical conductivity 
in the 2D orbital liquid regime (T > T*) can be compared with the conductivity in the 
3D-orbital liquid state. In particular there is (1) a broad continuum which decreases 
towards high frequency with a half-width uJi/2 ~ 2.5 and 3.5t in two and three dimen- 
sions, respectively. (2) The absorption increases with decreasing temperature and (3) 
the absolute value of the incoherent part of o"(co') at low frequency and temperature 
(Tq"''^ ~ 0.4 • 10^(ficm)~^ has the correct order of magnitude as in Tokura's experiments. 
(4) Despite this strong incoherence there is a finite Drude peak of about 20% of the 
total sum rule for small J. (5) The coherent motion and the Drude peak become more 
pronounced for larger values of J due to the 3-site hopping processes in the model. For 
the value J/t = 0.25 estimated here the relative Drude weight is 40 (20) % in the model 
with (without) 3-site hopping processes, respectively. (6) The low temperature values 



for the DC-conductivity and the resistivity can be estimated by assuming a width F 
for the Drude peak (taken from experiment). The additional scattering processes con- 
tributing to r are due to impurity or grain boundary scattering, i.e. these are extrinsic 
scattering processes which are not contained in the orbital model. The important infor- 
mation which follows from the model is the weight of the Drude peak, which together 
with r determines the DC-conductivity. 

The relatively small Drude peak in comparison with the incoherent part of cr(a;) 
clearly indicates that the carrier motion is essentially incoherent. Nevertheless we expect 
a quasiparticle peak in the single particle Green's function, yet with small spectral 
weight. Our study of the frequency dependent conductivity shows that the (saturated) 
ferromagnetic state in the manganites has unconventional transport properties due to 
the scattering from orbital excitations in combination with the exclusion of double 
occupancy. A detailed study of the orbital dynamics in the doped system is necessary 
for a deeper understanding of these issues. 

The orbital model resembles the t-J model, which describes the correlated motion 
of charge carriers in the cuprate superconductors. Yet an important difference is the 
cubic symmetry of the pseudospin representation of the orbital degrees of freedom. The 
hopping matrix elements depend on the orbital orientation, and are in general different 
in the two diagonal hopping channels. Moreover there is also an off-diagonal hopping 
matrix element, which does not exist at all in the t-J model. At low temperatures 
(T < T* 0.2t) doping induces (ferromagnetic) x^-y^ orbital order for realistic values 
for the exchange interaction J in the 2D orbital model. This is in striking contrast to 
t-J physics in cuprates, where the system changes upon doping from antiferromagnetic 
long-range order to an antiferromagnetic spin liquid. In both models coherent motion 
coexists with strong incoherent features. 

Our calculations have shown that the orbital mechanism can explain the order of 
magnitude of the conductivity at low temperature. The orbital degrees of freedom are 
certainly also important for a quantitative calculation of the coUossal magnetoresis- 
tance. This is obvious because in the full Kondo lattice model there is a close interplay 
between orbital and spin degrees of freedom. 

We acknowledge helpful discussions with F. Assaad, J. van den Brink, L. Hedin, G. 
KhaliuUin, A. Muramatsu, A. M. Oles and R. Zeyher. 
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